Nucleosynthesis in massive stars revisited 
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We have performed the first calculations to follow the evolution of all stable nuclei and 
their radioactive progenitors in a finely-zoned stellar model computed from the onset of 
central hydrogen burning through explosion as a Type II supernova. Calculations were 
done for 15 M , 20 M , and 25 M Pop I stars using the most recently available set of 
experimental and theoretical nuclear data, revised opacity tables, and taking into account 
mass loss due to stellar winds. Here results are presented for one 15 M model. 

1. INTRODUCTION 

Stars above ~ 10 M are responsible for producing most of the oxygen and heavier 
elements found in nature. Numerous studies have been devoted to the evolution of such 
stars and their nucleosynthetic yields, e.g., |l]||. However, our knowledge of both the 
input data and the physical processes affecting the evolution of these stars has improved 
dramatically in recent years. Thus, it became worthwhile to improve on and considerably 
extend the previous investigations of pre- and post-collapse evolution and nucleosynthesis. 
We present the first calculation to determine, self-consistently, the complete synthesis of 
all stable nuclides in any model for a massive star. Due to the limited space, in this 
report we mainly focus on giving an outline of our investigations. For further details on 
the calculations and results the reader is referred to the full papers ||[|] . 

2. REACTION NETWORK 

We employed a nuclear reaction network of unprecedented size in full stellar evolution 
calculations. The network used by |L] (WW95), large in its day, was limited to 200 nuc- 
lides and extended only up to germanium. Studies using reaction networks of over 5000 
nuclides have been carried out for single zones or regions of stars, especially to obtain 
the r-process, e.g., [H-H], but "kilo-nuclide" studies of nucleosynthesis in complete stellar 
models (typically of 1000 zones each) have been hitherto lacking. Similar to WW95, 
nucleosynthesis was followed by co-processing the stellar model throughout its evolution 
using the extended nuclear reaction network. From hydrogen ignition through central 
helium depletion a 617 nuclide network was employed that included all elements up to 
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Figure 1. Production factors in the ejecta of a 15 M star relative to solar abundance 
vs. nuclear mass number. 



polonium, adequate to follow the s-process. Just prior to central carbon ignition, we 
switched to a network containing 1482 nuclides. That network incorporated more neutron- 
rich isotopes to follow the high neutron fluxes in carbon (shell) burning and was also 
slightly extended on the proton- rich side to follow the 7-process The nucleosynthesis 
during the supernova explosion itself was followed in each zone using a 2437 nuclide 
network including additional proton-rich isotopes to better follow the 7-process in the 
neon-oxygen core, and also many additional neutron-rich isotopes to follow the n-process 
expected during supernova shock front passage through the helium shell. Here we will 
ignore the nucleosynthesis that occurs in the neutrino wind which may be the principal site 
of the r-process because its thermodynamic properties continue to be poorly understood. 



3. INPUT PHYSICS 

Our calculations were performed using the stellar evolution code KEPLER |l] with 
several modifications relative to WW95 (mass loss due to stellar winds ||10|| , improved 
adaptive network) and updates (OPAL95 opacity tables [PHI , neutrino loss rates ||T2f ). 
We generated a new library of experimental and theoretical reaction rates. As the basis 
of our reaction rate set we used statistical model calculations obtained with the NON- 
SMOKER code fll||l~4|| . A library of theoretical reaction rates calculated with this code 
and fitted to an analytical function — ready to be incorporated into stellar model codes 
- was published recently [f[4| . It includes rates for all possible targets from proton to 
neutron dripline and between Ne and Bi, thus being the most extensive published li- 
brary of theoretical reaction rates to date. For the network described here we utilized 
the rates based on the FRDM set. This was supplemented with experimental neutron 
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capture rates along the line of stability ||15|| . Experimental (0,7) rates were implemented 



for 70 Ge [|16j and 144 Sm [T^]. The derived a+ 70 Ge and «+ 144 Sm potentials were also 
utilized to recalculate the transfer reactions involving these potentials. For the important 
rates for a-capture reactions on self-conjugated (N = Z) nuclides, a new semi-empirical 
rate determination was implemented fT8|. For comparison, we used different sets of ex- 
perimental and theoretical rates for elements below neon: WW95, Ref. [|T^], and NACRE 
p0|| . Experimental (3~, f3 + , and a-decay rates were taken from []2T| and theoretical (3~ 
and /3 + rates from p2|. As a special case, we implemented a temperature-dependent 180 Ta 



decay p3| . For A < 40 we also included recent theoretical weak rates p4| . We did not 



follow the ^-process for nuclides with Z or iV larger than 40. The supernova explosion 
was simulated, as in [[!]], by an inward-outward moving piston resulting in a total kinetic 
energy of the ejecta at infinity of 1.2 xlO 51 erg. The final mass cut outside the piston is 
determined self-consistently from the hydrodynamic calculation. 

4. RESULTS AND DISCUSSION 

Here we only summarize the important results for one 15 M Pop I star. The production 
factors of this non-rotating model are shown in Fig. [I] as an example of our results. The 
model included mass loss and the rate set of [BJ below Ne. Note, however, that though this 
mass is a numerically typical case of a Type II or lb supernova, the average nucleosynthetic 
yield of massive stars is the result of populations of different stars each of which has its 
own peculiar yields which must be combined to result in a solar-like abundance pattern. 
Other stars and rate sets will be discussed in ||. 

4.1. Stellar structure 

The revision of the opacity table and the introduction of mass loss generally leads to 
smaller helium core sizes which tend to also decrease the mass of the carbon-oxygen and 
the silicon core. Note, however, that the absolute values of these core masses depend 
on the uncertainties, in particular, of the mixing processes in the stellar interior, such as 
semiconvection, overshooting, and rotation. 



The change in the weak rates |[24|| , important after central oxygen burning, leads to a 
2 — 3 % higher electron fraction per nucleon, Y e , at the time of core collapse in the center of 
the star and the "deleptonized core" tends to comprise less mass [f|]. More important for 
the core collapse supernova mechanism might be the 30 — 50 % higher densities of the new 
models between the region of m = 1.5 — 2 M© Q, which may result in a correspondingly 
higher ram-pressure of the infalling matter. 

4.2. Intermediate and heavy element nucleosynthesis 

A strong secondary s-process contribution appears between iron and a mass number of 
A = 90. Above A = 100 the s-process in our 15 M Q star is very weak, but it becomes 
notably stronger in stars with more massive helium cores that perform helium burning at 
higher entropies. Furthermore, the strength of the s-process is found to be very sensitive to 
the (a,n)-(a,7) rates and branching on 22 Ne which is experimentally not well determined 
(see also p5|). Second only to the well-known strong dependence of the stellar structure 
on the 12 C(a,7) rate, it becomes another important candidate for further laboratory study. 

The proton-rich heavy isotopes above A = 123 can be well produced by the 7-process 



4 



occurring during implosive and explosive oxygen and neon burning. The proton-rich 
isotopes around A = 160 and those between A = 100 and A = 123, however, are under- 
produced by a factor of 3 to 4 with respect to 16 O. The isotope 180 Ta seems to show a 
strong overproduction by the 7-process. However, in the figure only the totally produced 
180 Ta is shown. The surviving yield can be brought down to a more consistent production 
level by accounting for the distribution between ground and isomeric states M. 

The expected r- or n-process production by the supernova shock front passing through 
the base of the helium shell is not significant in any of our model stars, not even at 
A = 130. We observed some redistribution of isotopes at the base of the helium shell 
around A = 123 but this did not show the characteristics of a typical r-process nor was 
it important compared to the total yield of the star. 

Summarizing, we have presented the first calculation to follow the complete s-process 
through all phases of stellar evolution and the 7-process in the whole star through the 
presupernova stage and subsequent supernova explosion. This research was supported, 
in part, by DOE (W-7405-ENG-48), NSF (AST 97-31569, INT-9726315), the Alexander 
von Humboldt Foundation (FLF-1065004), and the Swiss NSF (2124-055832.98). 
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